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I— i! ABSTRACT 

^ . 

Context. The association of GRB 980425 with the nearby supernova SN 1998bw at z=0.0085 implies the existence of a 
^ population of gamma-ray bursts with an isotropic-equivalent luminosity which is about 10'* times smaller than in the standard 
^ cosmological case. Apart from its weak luminosity, GRB 980425 appears as a "normal" burst regarding all its other properties 
C/2 , (variability, duration, spectrum), with however a rather low peak energy Ep ~ 30 — 100 keV. 

/ Aims. We investigate two scenarios to explain a weak gamma-ray burst such as GRB 980425 : a normal (intrinsically bright) 

gamma-ray burst seen off-axis or an intrinsically weak gamma-ray burst seen on-axis. 
, Methods. For each of these two scenarios, we first derive the conditions to produce a GRB 980425-like event and we then discuss 
^ ' the consequences for the event rate. In the second scenario, this study is done in the framework of the internal shock model. 
T-H I Results. If we exclude the possibility that GRB 980425 is an occurence of an extremely rare event observed by chance during 

■ the first eight years of the "afterglow era" , the first scenario implies that (i) the local rate of "standard" bright gamma-ray 
' bursts is much higher than what is usually expected; (ii) the typical opening angle in gamma-ray bursts ejecta is much narrower 

than what is derived from observations of a break in the afterglow lightcurve. In addition to this statistical problem, we show 
f — ■ that the afterglow of GRB 980425 in this scenario should have been very bright and easily detected. For these reasons the 
' second scenario appears more realistic. We show that the parameter space of the internal shock model indeed allows GRB 
C , 980425-like events, in cases where the outflow is only mildly-relativistic and mildly-energetic. The rate of such weak events 

■ in the Universe has to be much higher than the rate of "standard" bright gamma-ray bursts to allow the discovery of GRB 
^ ' 980425 during a short period of a few years. However it is still compatible with the observations as the intrinsic weakness 

• i-H , of these GRB 980425-like bursts does not allow detection at cosmological redshift with present gamma-ray instruments. We 
' finally briefly discuss the consequences of such a high local rate of GRB 980425-like events for the future prospects of detecting 
^ , non-electromagnetic radiation, especially gravitational waves. 

Key words. Gamma rays: bursts - (Star:) Supernovae: individual: SN1998bw - Shock waves - Radiation mechanisms: non- 
thermal 



1. Introduction 

Our knowledge of the gamma-ray burst (hereafter GRB) 
population has dramatically increased in the recent 
years. In addition to "standard" cosmological bright 
GRBs peaking at 100 keV - 1 MeV have been dis- 
coved X-ray Flashes (hereafter XRFs) and X-ray Rich 
GRB s (hereafter XRRs) that peak at much lower en- 
ergy (jHeise et al.ll200ll: iKippen et 311120011 iBarraud et al.l 
20031) weak GRBs (such a s GRB 031203, ^azonOT_et al. 
2004ISoderberg et al.ll2004l) and even a few "local" events: 



GRB 98042 5 taala ma et al 
( Mirabal et al. 2006) 



1998h and GRB 060218 



It is usually believed that a long 
GRB is produced by an ultra-relativistic outflow ejected 
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from a black hole that has just form ed in the c ollapse 
of a very massive star (coUapsars, Wo oslevI (|l993 )). The 
prompt GRB emission originates from the relativistic 
ejecta itsel f, and is prob ably due to the formation of inter- 
nal shocks (|Rees fc Meszaros 1994) , whereas the afterglow 
is emitted by the strong shock that propagates within the 
ambient medium d uring the deceleration phase of the out- 
flow (Me szaros &: R ccs 1992). 

It is important to understand the physical origin of the ob- 
served diversity in this theoretical framework. There are 
two groups of explanations : either the diversity is only 
apparent: all GRBs are intrinsically very similar, but the 
observation conditions (viewing angle for instance) break 
this similarity; or the diversity is intrinsic: all coUapsars 
do not produce the same relativistic outflow and therefore 
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the same GRB. 

In this paper, we focus on the closest GRB, i.e. GRB 
980425. This burst has been detected by Beppo-SAX and 
BATSE. It appears as a "standard" single-pulse burst re- 
garding its gamma-ray properties : duration ^ 31 s, 
peak flux (40-700 KeV) P ~ 2.4 x lO"'^ erg/cmVs, Band 
spectrum with low- and high-energy slopes a —0.8 
and ^ ~2-3, peak energy Ep ^ 68 ± 40 keV 



Source at z 



( Frontera et al. 200G| ). With maybe the exception of a 
somewhat lower peak energy, all these properties are 
very close to thos e of standard GRBs. Note also that 
Norris et all (|2000h derived an especially large time lag for 
this burst between BATSE channels 1 and 3 : A13 ~ 4.5 s. 
A peculiar type Ic supernova, SN 1998bw, was dis- 
covered simultaneou sly in the Beppo-SAX error-box 
( Galama et al.l [l998h . The probability to have two such 
rare events occuring at the same time in the same direc- 
tion is very low, so that GRB 980425 and SN 1998bw are 



proba bly physically associated (see also lKouveliotou et al 



l2004l) . This is reinforced by other associations that have 
been found since 1998, the better case being probably the 



association of GRB 030329 with SN2003dh ( Stanek et al 



20031 ) at z = 0.168. If the association GRB 980425 / SN 
1998bw is real, then this burst is far from "standard" . The 
host galaxy of SN 1998bw is indeed located at z = 0.0085 
and GRB 980425 is therefore the closest GRB ever de- 
tected. As the peak flux of GRB 980425 is comparable to 
typical peak fluxes of other Beppo-SAX GRBs, it means 
that it is intrinsically much weaker (by more than four 
orders of magnitude). 

In section [3J we recall the various physical factors enter- 
ing in the observed peak flux and peak energy of a cos- 
mic GRB. We then study in section [3] a scenario where 
GRB 980425 is a standard bright GRB seen off-axis, and 
therefore apparently weak. In section [¥] we detail an al- 
ternative explanation, where GRB 980425 is intrinsically 
weak. We show how the internal shock model can allow 
for such weak bursts. Both in section [3] and 2] we discuss 
the consequences of each scenario in terms of GRB rate. 
Our results are summarized in section [51 

2. Peak flux and peak energy of a cosmic GRB 

We assume that a GRB is produced by a relativistic out- 
flow of Lorentz factor F and opening angle A9 ^ 1/F 
generated by a source at redshift z (see Figure [T|). We de- 
fine ^0 as the angle between the line-of-sight and the axis 
of the ejecta. The observed bolometric peak flux and peak 
energy are given by 



= /Cp (F; AO; Oq) 



and 



1^ 



En 



l + z ' 



(1) 



(2) 



where Lrad,4-n- and Ep are the isotropic equivalent lumi- 
nosity and the peak energy measured in the GRB source 




^e>l/^ 



~~. Obs 

Relativistic wind X " (on axis) 



■ Obs 
(off-axis) 

Fig. 1. Geometry. A relativistic outflow of Lorentz fac- 
tor F is emitted by the source (located at redshift z) in a 
cone of opening angle A6'. The line-of-sight of the observer 
makes an angle 9q with the axis of the cone. The obser- 
vation is made on-axis (resp. off-axis) if ^ A0 (resp. 
6*0 > Ad). 

frame. The correction for the viewing angle are approxi- 
matively given in the on-axis and the off-axis cases by 




1 if 6I0 < AO 

-, i if 6I0 > AQ 

2(i+r2(0o-Ae)^)-' 



/Ce (F;A0;0o) = 



i+r2( 



if 6*0 < AB 
if 00 > A^? 



(3) 



(4) 



Therefore, the peculiar properties of GRB 980425 may 
have two origins : (i) either GRB 980425 is a "stan- 
dard" GRB, i.e. intrinsically bright with irad.47r ^ 
a few 10^'^ erg.s"^ , but appears as a weak GRB because 
it is seen with a large viewing angle; (ii) or GRB 980425 
is an intrinsically weak GRB (Lrad,4ir — 3 x 10'*® erg.s""'^) 
seen on-axis. In this case GRB 980425 has been detected 
only because of its very low redshift. In the following, 
we successively investigate these two scenarios. The first 
one has already been adressed by several authors (e.g . 
Nakamuralll999l: ISalmonsonll200ll lYamazaki et al.ll20oi 
Guetta et al.l 2004[ l so we focus mainly on the second sce- 
nario. For this case, the intrinsic properties of GRBs are 
studied in the framework of the internal shock model, 
where the prompt emission comes from a relativistic wind 
which converts part of its kinetic energy into radiation via 
the formation of shock waves within the wind itself. Such 
internal shocks can occur if the wind is generated with a 
highly variable Lorentz factor (Rees fc MeszaroSiil994) . 

3. GRB 980425 as an intrinsically bright GRB 
seen off-axis 

3.1. Conjugate effect of a low redshift and a large 
viewing angle 

In Figure[2l GRB 980425 has been plotted in a bolometric 
peak flux vs peak energy diagram, together with the 
GRBs detected by Beppo-SAX which have a known 
redshift. A sensitivity of 10~^ erg cm~^ s~^ represen- 
tative of the BATSE and Beppo-SAX GRBM+WFC 
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Fig. 2. Bolometric peak flux vs peak energy di- 
agram : All Beppo-SAX GRBs with a known redshift 
are indicated in this diagram (crosses). The peculiar 
burst GRB 980425 is indicated by a circled cross. A 
vertical dotted line stands for a threshold Pmin = 
10~^ erg cm"^ s^^ representative of the BATSE an d 
Beppo-SAX WFC+GRBM instruments |Bandl |2003| ). 
Effect of the redshift : we have plotted the evolution of 
GRB 980425 when its redshift increases from z = 0.008 
to z = 1 (cosmological distance). The final result is indi- 
cated as GRB 980425*. Effect of the viewing angle : we 
have then plotted the evolution of GRB 980425* when the 
viewing angle decreases. We assume that GRB 980425* is 
seen off-axis with Oq = A0-|-4/r (see text). The final result 
for 6*0 = (on-axis) is indicated as GRB 980425**. 



instru ments is also indicated as a dotted line (jBandl 
20031 ). We use equations ©-(HI) to estimate under which 
conditions on the redshift and the viewing angle GRB 
980425 can be an intrinsically bright GRB despite its 
low apparent intensity. To clarify the respective roles of 
redshift and viewing angle, we perform the following two 
transformations : 

- Effect of the redshift : GRB 980425 has first been 
moved from z = 0.008 to z = 1 (cosmological distance). 
Both the peak energy and the peak flux decrease due 
to the increase of redshift and luminosity distance. The 
final result is named GRB 980425* : this burst is clearly 
much to weak to be observed as a GRB or even a XRF. 
Its bolometric peak fiux is indeed at least three orders 
of magnitude below the sensitivity of past or current 
detectors. 

- Effect of the viewing angle : the viewing angle of 
GRB 980425* is then decreased down to = (on-axis 
observation), assuming that GRB 980425 is initially seen 
off-axis with a large viewing angle = + k/T (here 



Fig. 3. The ratio of detected ofF-axis to on-axis 
GRBs in the local Universe : this ratio within a dis- 
tance D is plotted as a function of D for a power-law lumi- 
nosity function of slope —1.7 between Lmin = 10^" erg s~^ 
and Liiiax = 10^** erg s~^ and three different distributions 
of the opening angle: (1) a uniformly distributed open- 
ing angle between 1 and 10° (solid line); (2) a power-law 
distribution p(A0) cx A0^^ between l/F and 7r/2 (dashed 
line); (3) an opening angle correlated with the isotropic 
luminosity (dotted line). In this last case, the true lumi- 
nosity L{1~ cos AO) is assumed to be constant and equal 

to £min- 



fc = 4). This evolution is computed from the approximate 
formulae ^ and (H]). The final result is named GRB 
980425**. Now this burst is clearly back in the "standard" 
GRB region, very close to GRB 971214. 
Therefore, as it was already pointed out by 

the 



Yamazaki et al.l (|2003l ). the peculiar properties of 
GRB 980425 are compatible with those of an intrinsically 
bright GRB observed off-axis. In the case illustrated in 
Figure m the decrease of the flux of GRB 980425 by 
a factor ~ 10^ due to a viewing angle 6*0 = -f 4/r 
is compensated by a smaller luminosity distance 
{Di,{z = Qm8)/Di^(z = l)f - 2 X lO-'^. 

If it is clear from the previous analysis that a slightly 
off-axis viewing angle [Oq - A9 ^ 4/r - 2.3° if F ~ 100) 
can account for the properties of GRB 980425, one should 
however keep in mind that the peak flux decreases very 
rapidly with viewing angle (see equation ([3])). For 0q — 
A9 + fc/F with fc > 6, the peak flux is divided by a factor 
larger than ^ 10*^ and the decrease of luminosity distance 
is not able anymore to compensate : such off-axi s GRBs 
cannot be detected. As F ^ 100 in GRB outflows ijBaring 
19951: iLithwick k Sari|[200lh . the angle 6/F ^ 3.4° is smafl 
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and the probability to observe GRB 980425-like events 
should be low. 



3.2. A statistical problem ? 

For a detector having a threshold Pmin (hereafter we as- 
sume Pmin = 10~^ erg cm~2 s~^, which is representa- 
tive of instruments such as Beppo-SAX GRBM+WFC 
or BATSE) the observed rate of on-axis GRBs within 
D = AO Mpc is given by 



47r 

^on = -irD^o X / dA6' p (AO) (1 - cos A9) 
3 Jo 



(5) 



where Rq is the local GRB volumic rate and p{A9) is the 
probability distribution of the opening angle A9. Here we 
assume that the minimum GRB luminosity Lmin is larger 
than 47rl?2p,nin = 2 x 10^^ erg s~^ so that all nearby on- 
axis GRBs are detected. On the other hand, the rate of 
off-axis GRBs in the same volume is affected by the rapid 
decrease of the GRB flux with the opening angle so that 



TZoS = Ro dD 4ttD^ 



■it/2 



X COS A0 - cos A0 



dA6ip(A6i) 
.x(Ai) 



dL p{L) 
(6) 



where p{L) is the GRB luminosity function and fcmax is 
given by 



fcmax (-D, L) 




2 AnD^Prr 



1/3 



1 



(7) 



In both equations ^ and ([7]) , the probability distribution 
of the Lorentz factor F is neglected and F is supposed 
to be constant. In the most simple case where both the 
luminosity L and the opening angle A0 are constant, we 
find that the ratio of detected off-axis over on-axis GRBs 
can be expressed as a ratio of two solid angles : 



V 



COS A9 — COS A^ 



1 - COS A0 



(8) 



For L) = 40 Mpc and L = 10''^ erg s~\ we get fcmax = 5.4, 
so that fcmax/F <, 3.1° for F ^ 100. This leads to a rate 
of detected off-axis over on-axis GRBs TZos/TZon ^ 16 for 
A0 = 1°, 7^off■/7^o„ 1 for A0 = 7.5°, and 7^off/7^on 
0.7 for A0 = 10°. This implies that the rates of detected 
off-axis and on-axis GRBs should be comparable for typ- 
ical A0 of a few degrees and that it is only for A0 ^ 1° 
that the number of detected off-axis GRBs is much larger 
than the number of on-axis GRBs. The local apparent 
rate of "standard" on-axis bright GRBs can be estimated 
to be of the order of ^ 1/(5 000 - 30000 y r) within 40 
Mpc (IPorciani fc Madaull200ll[Schmidt 200ll : IPerna et al 



120031 : f Guetta et al.ll200^ . This simple analysis would then 



imply that with GRB 980425 we have observed by chance 
a very rare event. Even if it has some rather different prop- 
erties, GRB 060218/SN2006aj at z = 0.0331 (144 Mpc) 



( Cusumano et al.ll2006l: iMirabal et al.ll2006l : [Masetti et aP 
2006 : ISoderberg et al.j2006a )) is another case of a nearby 



burst. At this distanc e 7?.off/7?.nn shoul d be even smaller. 



Statistical studies (Bosniak et al. 20061 ) also indicate that 
there are probably other GRB980425-like events in the 
BATSE catalog. For all these reasons, we conclude that 
the off-axis interpretation probably suffers a statistical 
problem. In the more general case, the two solid angles 
in equation ([5]) have to be averaged over the luminos- 
ity function, the distance and the opening angle distri- 
bution, according to equations ([6]) and ([7]). We did that 
for a power-law luminosity function with slope —1.7 (this 
slo pe giving a good fit to the log N — log P diagram, see 



e.g Daigne et al. 2006f ) between Lmin — 10^° erg s ^ and 



Lmax = 10 erg s~ . We tested three possible distribu- 
tions for the opening angle: (1) a uniformly distributed 
opening angle; (2) a power-law distributed opening angle; 
(3) an opening angle correlated with the isotropic equiva- 
lent luminosity so that the true luminosity is constant, as 
suggested by obs ervations of achr omatic breaks in after- 
glow lightcurves ([Frail et al. I l200l[ ). We find (see Figure [3j) 
that except below ^ 5 Mpc where the total (on-|-off axis) 
event rate is very low due to a small volume, the expected 
rate of off-axis GRBs is either comparable (uniformly dis- 
tributed opening angle) or smaller than the on-axis GRB 
rate. This is in full agreement with our simple estimate 
(equation ([8])). This equation also indicates that the only 
way to have a much higher rate of off-axis GRBs is to 
assume a very small opening angle (AO <^ l /F), in con- 



tradiction with observations ([Frail et al.ll2001[ ) 



4. GRB 980425 as an intrinsically weak GRB 

We now consider in this section an alternative scenario 
where GRB 980425 is an intrinsically weak GRB seen on- 
axis. The intrinsic GRB properties are considered in the 
framework of the internal shock model ( Rees fc Meszaro^ 



4.1. Internal shocks 

The dynamics of the internal shocks as well as the tem- 
poral and spectral pr operties of the emission have been 
studied in details in Daigne fc MochkovitchI ( 1998h us- 
ing a simplified model where the outflow is made of 
a large number of discrete relativistic shells that in- 
teract by direct collision only. This approach has been 
validated by ID relativistic h ydrodynamical simulations 
( Daigne fc Mochkovitch 2000f ). We adopt here an even 
simpler version of the model, where we consider only a 
typical internal shock due to the collision of two shells of 
equal mass M and Lore ntz factors Fi and F2 > Fi. This 
toy model is described in lBarraud et al.l ( 20051 ). It is found 
that the collision between the two shells occurs at radius 



Rk 



(9) 
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that the isotropic equivalent radiated luminosity due to 
internal shocks is given by 



K + 1 

and that the peak energy takes the form 



^2x p6a:-l 



(10) 



(11) 



where E is the isotropic equivalent kinetic energy injection 
rate in the relativistic outflow, r is the duration of the 
relativistic ejection (i.e. here the time interval between the 
two shell ejections), F is the mean Lorentz factor (Fi + 
F2)/2 and the contrast k = r2/Fi is a measure of the 
initial amplitude of the variations of the Lorentz factor in 
the ejecta; £c is the fraction of the dissipated energy in the 
shock that is injected into non-thermal electrons, and K 
is a constant. The function (pxyi^) is given by 



4>xy{n) 



(12) 



and is steadily increasing with k. The x and y parameters 
have been introduced by considering that the peak energy 
in the comoving frame of the shocked material scales as 



E'p cx p^e^ 



p is the comoving density and e the dissi 
energy per unit mass. 



where 



lissipated 



In the standard synchrotron model with constant 
equipartition parameters it is assumed that the mag- 
netic field is amplified to reach a fraction ee of the 
total dissipated energy in the shock, and that a frac- 
tion 60 of this dissipated energy is injected into a frac- 
tion C of the electrons, which arc therefore accelerated 
to high Lorentz factors. This leads to x — 1/2 and 
y = 5/2, as the synchrotron energy scales as E'^ cx BT^, 

the magnetic field as B cx e^^{pe)^/'^ and the typi- 
cal Lore ntz factor o f the electr ons a s Fe oc {ec/Qe. 
However Daigne fc MoclikovitchI (|2003h have shown that 
smaller values of x and y are required to reproduce 
the hardness-intensity and hard ness-fluence correlations 
observed in many GR B pulses (jGolenetskii et al. 
Liang &: Kargati 3 ll99fih . They consider x = y = 0.5 and 
X = y = 0.25. Such exponents can for example be obtained 
with the standard synchrotron process if the equipartition 
parameters vary with shock intensity. 



with T = 5 s (observed duration of (1 + z)t = 10 s), 

= 10^1 erg.s^i) 



= 4, £; = 1.5 X 10 



^2 erg.s ^ {Ly 



ad,4ir 



and F = 300 has an observed peak energy E° — 200 keV. 
Note that in this case, the efficiency of internal shocks is 
Li-ad,4Tr/-£' — 6.7%, as the product of Cc — 1/3 with an 
internal shock dynamical efficiency of 20% for k = 4 (see 
eauation llOl). The other parameters are li mited by several 
constraints ( Daigne fc MochkovitchI 2004h : 



(1) Transparency during the internal shock phase. The 

relativistic ejecta has to be transparent during the internal 
shock pha se. The Thomson optical depth of the outflow is 
given by i Meszaros fc Rees 2000l : Daigne fc MochkovitchI 
2002, 1 : 



T(i?is) 



gtEt 



47rFi?is (i?is + 2F2 ct) m^c^ 



Then the condition T{Ris) < 1 leads to 



E 47rTOnC* 



F5 
with 



< 



/i(k) r , 



(13) /i(«) 



(k-1)(k+1)3 V(k-1)(k + 1)^ 



(14) 



(15) 



(16) 



(2) No pair production. The relativistic ejecta has to 
be transparent to pairs during the internal shock phase. 
Pairs are produced by photon-photon annihilation in 
the high-energy part of the GRB spectrum. At radius 
Ris, the approxirnate optical depth for pair creation is 
(|Meszaros fc Reesll2000tTLithwick fc Sarill200ll) : 



T±{Ris) 



a± 



C^TLradATrT ( Ep 



A'KR'^TmcC^ yFmcC^ 



f3-2 



(17) 



where [3 is the high-energy slope of the spectrum and a± 
is a dimensionless number that depends on the spectral 
shape. In the following, we adopt fi = 2.5 (|Preece et al. 
,2000). This leads to a± ~ 0.06. Then the condition 
T±{Ris) < 1 leads to 



Y5+6x(J3-2) 

with 



a±aT 



K 



2-/3 



,l+2x(/3-2) 



(18) 



4.2. Producing GRB980425-like events 

The value of I/rad,47r and Ep given by equations (fTO|) and 
(fTTj) are fixed by 6 physical quantities : E, F, k, r, Cc and 
K. The value of eo and K depend on the details of the 
physical processes in the shocked material, which are not 
studied here. We adopt = 1/3 and we fix the value of 
K by demanding that a "typical" GRB at redshift z = 1 



/2(«) 



64k4 



(y^-l)2(K-l)2(K-f 1)5 $^;2(k) 



(19) 



We also considered an additional constraint provided 
by the lack of an observed high-energy cutoff due to 
photon-photon annihilation in GRB spectra (see e.g. 
iLithwick fc Sari|[2001,') . However, the constraint expressed 
by equation ([TS]) is dominant in most cases; 
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Fig. 5. The internal shock model parameter space 

at z = 1. The diagram plotted in Figure [4] has been 
moved from z = 0.0085 to z = 1, k, r, Cc and K be- 
ing kept constant. Only the case where a; = y = 1/4 and 
the external medium has a constant density n — 1 cm~'^ 
is shown. Weak GRBs are now undetected. The thick line 
corresponds to Ep = 150 keV. Above this line are found 
Ep — 400 keV and 1 MeV. Below are plotted lines for 
Ep = 70 keV, 40 keV and 10 keV. The big dot indicates 
the location of a "standard" GRB with Ep = 150 keV 
(observer frame) and Lrad,47r — 10^^ erg s^^. 



(3) Internal shocks operate before the deceleration 
phase. The internal shock phase has to occur before the 
deceleration of the cjccta due to the external medium. For 
a density p oc A/R^ in the environment, the deceleration 
starts at 



^dec 



Et 



Then the condition i?is < i?dcc leads to 



E 47r ^ 
> Ad 



(k- 1)(k + 1)S 



3-s 



(20) 



(21) 



Figure [3] shows the location of these constraints as 
well as lines of constant peak energy in the T-E plane 
for two sets of the x and y parameters : x = y — 1/4 
(top) and X = y = 1/2 (bottom) and for two possible 
environments : either a uniform medium (s = 0) with 
n = A/nip = 1 cm~^ (left column) or a dense stellar 
wind (s = 2) with A, = A/ (5 x 10" g cm^^) = 0.1 
(right column). The other parameters are fixed : r = 30 s 
which is the duration of GRB 980425 and k = A. The 
transparency constraint is dominated by equation ([T8| 



(optical depth for pair creation) and not by equation \W\ 
(Thomson opacity of ambient electrons) . The environment 
constraint (equation (HI])) is much stronger in the case of 
a stellar wind. For winds denser than A^, =0.1 (WR stars 
may have A, ^ 1), it becomes very difficult for internal 
shocks to operate before decele ration starts. This p r oblem 
was already mentioned in iDaigne fc Mochkovitehl (|l999L 
200lh . 



Despite these constraints, the internal shock model can 
explain a large range of bolometric luminosities and peak 
energies, in agreement with the observed GRB diversity. 
As can be seen in FigurelU the internal shock model allows 
the existence of GRB980425-like events at z = 0.0085, 
except in the case of a wind environment with x = y = 
1/2. These weak GRBs correspond to mildly relativistic 
(r ~ 10 - 20) and mildly energetic (£; ~ 5 x lO'^'^ erg s^^) 
outflows. On the other hand Figure [5] confirms that at 
cosmological distance {z = 1), only "standard" GRBs with 
highly relativistic (F ^ 30) and highly energetic {E ^ 
1q5i gj.g g-i gg^j^ detected. Notice that Figure[5]show no 
X-ray Flashes (XRFs) or X-ray rich GRBs (XRRs), i.e. no 
GRB s with low peak ene rgies {Ep ^ 50 keV). As explained 
in lBarraud et al.| (2005), such soft GRBs are produced in 
outflows with high Lorentz factors F and low contrasts 
K, i.e. outflows with a very low baryonic pollution and a 
smoother ejection. This is illustrated in Figure [6] where is 
shown the effect of a low constrast k. 



4.3. The afterglow of GRB 980425-like events 

The afterglow of a "standard" GRB at 40 Mpc should 
be ~ 10 = 51og(D(z ^ l)/40 Mpc) magnitudes brighter 
than for a cosmological GRB, due to the short distance. It 
would then easily peak at the fifth magnitude in V band. 
Such a bright afterglow was not detected in association 
with GRB 980425. The emission in V ba nd in entirely 
domi nated by the lightcurve of SN1998bw ( Galama et al.l 
Il998 ) (see filled circles in the bottom panel of Fig. [7]) 
and no afterglow was detected. In X-ray, the WXM of 
Beppo-SAX has detected a weak decreasing source that 
could be the X-ray afterg low : see the fill ed circles in 
the top panel of Fig. [7] (jPian et al.l l2000l ). Note that 
only the first point is a firm detection. Any scenario for 
GRB 980425 should then predict a weak afterglow to be 
consistent with these data. We have plotted in Fig. [7] 
the predicted lightcurve in X-rays and in the V band of 
the afterglow for each of the two scenarios considered in 
this paper. Again, the external medium is either a stellar 
wind with s = 2 and A^ =0.1 (solid line) or a uniform 
medium with s — and n — 1 cm^'^ (dashed line). The 
first case (wind environment) should be preferred due 
to the association of GRB 980425 with SN1998bw. The 
afterglow is computed assuming "typical" values for the 
parameters : p = 2.5 (slope of the electron power law 
distribution), Cc = 0.1 (fraction of dissipated energy 
injected in relativistic electrons) and eb — 10^^ (fraction 
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Fig. 6. Effect of the contrast k. The diagram at z = 1 
plotted in Figure [5] has been moved to the situation where 
the amplitude of the initial variations of the Lorentz factor 
is low : K = 1.5 instead of k = 4. All other parameters are 
the same as in Figure [51 The thick line corresponds to 
Ep = 40 keV. Above this line are found Ep = 70 keV 
and 150 keV. Below is plotted the line for Ep — 10 keV. 
The big dot indicates the location of a "standard" XRR 
with Ep — 40 keV (observer frame) and irad,4ir = 3 x 
lO^^ergs-^ 

of dissipated energy injected in magnetic field). The 
thick lines correspond to an intrinsically weak GRB 
seen on-axis. According to our analysis of the prompt 
emission, we have adopted Tq = 10 for the initial Lorentz 
factor and Eq = 6 x 10'*^ erg for the initial kinetic energy. 
The thin lines correspond to an intrinsically bright GRB 
seen off-axis. The initial Lorentz factor and kinetic energy 
have "standard" values Fq = 100 and Eq = lO^^ ^j^^ 
opening angle is fixed to A9 = 5° and the viewing angle 
to 9o = A9 + 4/Fo (see sectional). 

It appears clearly that in the off-axis scenario, the pre- 
dicted afterglow is much too bright. It should have been 
easily detected in X-rays or in the V ba nd. This is again an 
evidence against the off-axis scenario (|Waxmanll2004allbl ). 
A discussion ba sed on the radio afterglow leads to th e 
same conclusion ( Waxmanll2004bl : Soderberg et al.| [2006aV 
To decrease the afterglow flux, one should either consider 
a very low density environment : A* ^ 3 x lO^"* (.s = 2) 
or n 10~^ cm~'^ (s — 0), or assume very low cquiparti- 
tion parameters such as eo < 5 x 10"'' or cb < 10^''. On 
the other hand, in the case of an intrinsically weak GRB 
seen on-axis, the afterglow is also very weak, due to the 
low value of the kinetic energy. The afterglow lightcurve 
in the V band is always well below the lightcurve of SN 



Fig. 7. Afterglow lightcurve in X-rays (top panel) 
and V band (bottom). The thick lines show the pre- 
dicted afterglow lightcurve in the case of an intrinsically 
weak GRB seen on-axis and the thin lines in the case 
of a bright GRB seen off-axis. The solid lines correspond 
to a stellar wind environment and the dashed lines to a 
uniform medium (see text). Observations are shown with 
filled circles. 

1998bw, while in X-rays, it is very close to the Beppo-SAX 
detections or upper limits. 

4.4. Rate of GRB 980425-like events 

The parameters of the internal shock scenario should 
in principle be given by models of the central engine. 
However, the production of the relativistic outflow is far 
from being understood. Therefore, the probability den- 
sity of parameters such as the injected kinetic energy 
flux E, the mean Lorentz F, the constrast k and the du- 
ration of the relativistic ejection r are unknown, which 
prevents us from using the internal shock model to pre- 
dict the rate of a specific class of GRBs, such as GRB 
980425-like events. However we have shown that in the 
framework of the internal shock model, "standard" GRBs 
are produced by highly energetic / highly relativistic out- 
flows, whereas GRB 980425-like events are produced by 
mildly energetic / mildly relativistic outflows. The fact 
that GRB 980425 has been observed (and maybe a some- 
how similar event with GRB 060218) whereas no "stan- 
dard" GRB has been observed within a few hundreds 
Mpc would then indicate that less spectacular relativis- 
tic ejections are more probable than the extreme outflows 
that are needed to produce a "standard" GRB. Indeed, 
the fraction of stellar collapses that produce bright GRBs 
is estimated to be of the order of 10"^ (apparent rate) 
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(iPorciani k Madaul [2OO1I ) . i.e ^ 10""^ 
rection for the b e aming angle ( Frail et al. 200l[) . Recently 
Soderberg et al. (|2006bl) have estimated that the rate of 
"sub-energetic GRBs" (GRB980425-like events) has to be 
^ 10 times larger than the rate of "standard" bright 
GRBs, which is still compatible with only a fraction of 
type Ic supernovae being associated with such GRBs if 
the beaming angle is not too small. 



4.5. Non electromagnetic emission 

As the rate of GRB 980425-like events is expected to 
be larger than the rate of "standard" GRBs, it is in- 
teresting to investigate whether such weak bursts could 
contribute in a significant way to the production of non 
electromagnetic radiation, i.e. ultra-high energy cosmic 
rays (UHECRs), high energy neutrinos and gravitational 
waves. Due to a lower Lorentz factor and a lower kinetic 
en ergy, it can be s hown, using the formalism developped 
bv lWaxmanI (|200lh that GRB 980425-like events are prob- 
ably not able to accelerate protons at energies higher than 
~ 10^^ eV. In addition, the low power of these bursts is 
not compensated by the larger event rate so that the local 
energy injection rate in UHECRs due to GRB 980425-like 
events is probably only a few percents of the rate due 
to "standard GRBs" . We conclude that GRB 980425-like 
events are not good candidates as sources for UHECRs, 
and probably do not even contribute significantly to the 
observed CR spectrum below 10^^ eV. The same kind of 
conclusion is reached for high-energy neutrinos. Gamma- 
ray photons of ~ 100 keV will interact with protons at 
the photo-mes on threshold of the A resonance (see e.g. 
Waxmanlliool . i.e. of energy 2 x 10^^ eV in the observer 



frame for F = 10. Due to the decay of the produced pions, 
neutrinos are emitted, which carry typically 5 % of the ini- 
tial proton energy, i.e. ^ 10^^ eV. However, the neutrino 
flux should a priori scale with the high energy proton flux, 
and therefore the contribution of GRB 980425-like events 
should be only a small fraction of the total GRB contri- 
bution. 

The prospects are much better for gravitational waves. 
Long GRBs are believed to be asso ciated with ma ssive 
stars that collapse into black holes (jWooslevI 119931 ). As 
only a fraction of these collapses lead to the production 
of a gamma-ray burst, peculiar conditions are required. It 
is natural to expect that a high rotation of the progenitor 
favors GRBs, as the ultra-relativistic outflow is probably 
easier to eject along the rotation axis. Therefore, coUap- 
sars might be asymetric and then be the source of gravita- 
tional waves. This is also supported by the physical mod- 
elling of SN19 98bw, that probably requires a high degre e 



emng ot &iNiiji j»Dw, tnat proDapiy requires a mgn aegre ( 
of asymmetrv (iHoflich et al.l 1999t iNakamura et al" 2001) 
iKobavashi &: Meszarod ( 2003 ) estimate the expected sig- 
nal for an asymmetric stellar collapse to be of the order 
of /i ~ 10-22 ^Yie 100-1000 Hz range at ~ 30 Mpc for 
plausible parameters. Such a signal will be detected by ad- 
vanced gravitational waves detectors. It is interesting to 



note that in the proposed scenario the nature of the initial 
event responsible for GRB 980425 may not be so different 
than for "standard" GRBs. This is supported by the sim- 
ilarities between SN2003dh (associated with GRB 030329 
at 2; = 0.168) and SN 1998bw. So the emission of gravita- 
tional waves of GRB 980425-like events may be compara- 
ble with those expected for "standard" long GRBs. Due 
to their higher rate within a few 10 Mpc, GRB 980425-like 
events are ther efore ve ry promising sources for advanced 
GW detectors (|Norri£ll200a) . 



5. Conclusion 

We have investigated two scenarios to explain GRB 
980425-like events : 

- in the first scenario, GRB 980425 is a normal (intrinsi- 
cally bright) GRB seen off-axis. A viewing angle that is 
greater than the opening angle of the jet by a few 1/F is 
enough to produce an apparently weak GRB like 980425. 
However this scenario suffers several problems : (1) the 
afterglow of GRB 980425 should have been very bright 
and therefore detected; (2) the rate of GRB 980425-like 
events should be very low. If we exclude the possibility 
that we observed by chance an extremely rare event and 
if we take into account the fact that a local on-axis intrin- 
sically bright GRB has never been observed, we get two 
strong constraints : the local rate of "standard" bright 
GRBs has to be much higher than what can be simply 
extrapolated from the cosmic event rate and at the same 
time the typical jet opening angle in GRBs must be much 
narrower than what is usually inferred from observations 
of breaks in afterglow lightcurves. 

- Since this appears rather unrealistic, we prefer the sec- 
ond scenario, where GRB 980425 is an intrinsically low- 
luminosity GRB seen on-axis. We have shown that the 
parameter space of the internal shock model allows such 
events, that would be produced by mildly relativistic 
(F ~ 10 - 20), mildly energetic (i; ~ 5 x 10"*^ erg.s"'') 
outflows. The consequences are that (1) in coUapsars, 
the central engine in most cases fails to power a highly- 
relativistic/highly-energetic ejection (necessary for bright 
"standard" GRBs) but can more easily power a mildly 
relativistic/mildly-energetic outflow; (2) the rate of GRB 
980425-like intrinsically weak events in the Universe is 
then much higher than the rate of "standard" intrinsically 
bright GRBs but the apparent rate is lower because these 
events can be detected only at low redshift. Finally, as 
this scenario infers a higher local rate for GRB 980425-like 
events compared to "standard" bright GRBs, we have dis- 
cussed the possible non-electromagnetic emission of these 
bursts. We find that they are not promising sources of 
UHECRs and high-energy neutrinos, as the correspond- 
ing flux probably scales with the gamma-ray flux. On 
the other hand, since GRB 980425 was associated with 
an extreme type Ic supernova as those associated with 
"standard" long GRBs, one can expect the corresponding 
stellar collapse to be asymmetric and GRB 980425-like 
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events could therefore be promising sources of gravita- 
tional waves within a few 10 Mpc. 
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Fig. 4. The internal shock model parameter space at z = 0.0085. In the T-E plane, lines of constant peak energy 
(observer frame) are plotted for k = 4, r = 30 s, ec = 1/3 and a fixed value of K computed as explained in the text. 
These lines are limited on the left side by the transparency constraint given by equation IjlSp (optical depth for pair 
creation : solid line) and ([TS]) (opacity due to ambient electrons : dashed line) and on the right side by the environment 
constraint given by equation (PT|) . As k and are fixed, the radiative efficiency is also fixed to 6.7% (see text) : thin 
dotted horizontal lines show the radiated luminosity L^ndATr corresponding to the injected kinetic energy rate E. The 
shaded area excludes GRBs that cannot be observed by a bolomctric detector of threshold Pmin ~ 10^^ erg cm^^ s^^. 
The thick line corresponds to Ep = 70 keV. Above this line are also plotted lines for Ep — 150 keV, 400 keV, 1 MeV, 
and 3 MeV. Below are found Ep = 40 keV and 10 keV. The location of GRB980425 is shown as a big dot. 



